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ABSTRACT 

New experimental results from photoproduction of hadrons at HERA are re- 
viewed. 

1. Introduction 

The photoproduction of hadrons was in the first years of HERA experimentation 
the process to be studied most easily because of the large cross section. The recent 
results tend to investigate the domain of harder physics which can be compared with 
QCD calculations. In this contribution the latest results are reviewed which were 
presentedB on the spring conferences of 1997 and are not covered by more specialized 
talks of this workshop. Specifically, I shall talk about inclusive particle production, 
resonances in multi-photon final state and prompt photons. 

2. Inclusive particle production 

The inclusive production of single hadrons in fixed target experiments and collid- 
ing beam experiments has been one of the important testing grounds for the QCD 
parton model. The HERA collider makes possible such studies in ep interactions 
over a wide range of squared four- momentum transfer, — Q 2 , from photoproduction 
Q 2 pa to very high photon virtuality. Recent progress in the determination of parton 
fragmentation functions enables us to make specific predictions for inclusive particle 
production in next-to-leading order QCD. The parton distributions of the photon and 
the proton constitute an important input to the cross section calculations. 



*To appear in the Proceedings of the Ringberg Workshop: New Trends in HERA Physics, Ringberg 
Castle, Germany, 25-30 May 1997. 



2.1. Neutral strange particles 

New results for K° and A inclusive photoproduction were obtained in the HI 
experiment! for an average jp cms energy W of 187 GeV. Kg mesons and A baryons 
are identified through their decay channels K s — > 7r + 7r~, A(A) — > p7i^(p7i + ), where 
the pion, the proton and the antiproton are reconstructed in the central jet chamber 
(CJC). In order to ensure an optimal acceptance of the CJC for neutral strange 
particles , their pseudorapidity has been restricted to \r)\ < 1.3 and their transverse 
momentum to 0.5 < p t < 5 GeV for K s and to 0.6 < p t < 5 GeV for A. All 
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Fig. 1. Measured K° (a,b) and A (c,d) total ep cross section in photoproduction as a function of pt 
and pseudorapidity r\ compared to PYTHIA in the kinematic range < 1.3 and 0.5 < pt < 5 GeV. 



results are multiplied by 2 and therefore correspond to K° + K° production and are 
denoted by K°. As the A cross section is within errors equal to A cross section they 
are combined and referenced as A. Fig. [I] shows a comparison of the measured K° 
and A cross sections in p t and rj with the predictions of PYTHIA in combination 
with JETSET using standard parameter settings (A s = 0.3 as well as parameter sets 
obtained from fits to DELPHI and E665 data)? In the case of K° mesons the measured 

a The following JETSET hadronization parameters have been changed with respect to the default 



2 



cross sections are in reasonable agreement with the Monte Carlo prediction using the 
DELPHI or E665 settings. However the measured A cross section lies significantly 
above the Monte Carlo prediction for all three settings. The UA5 collaboration has 
reported! a similar discrepancy between A production in pp interactions at 200 GeV 
and PYTHIA predictions. 

Recent NLO QCD calculations of inclusive spectra in ep scattering use fragmenta- 
tion functions! which have been fitted to e + e~ data. A comparison of this calculation 
to data in the region of p t > 1.8 GeV is satisfactory!'!. If fragmentation is process 
independent, the rate of low p t particles in jp and deep inelastic scattering (DIS) 
interactions will be the same. A comparison of K° production rates in 7p and Dis! is 
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Fig. 2. K° rate as a function of pt and centre of mass rapidity y* (upper part) and A rate as a 
function of y* (lower part) for non-diffractive events, jp data are represented by open and DIS data 
by closed circles. The inner vertical error bars indicate statistical error, the outer error bars the 
statistical and systematic error added in quadrature. 

shown in Fig. |2| (upper part) and of A production rates in Fig. || (lower part) for the 
non-diffractive event samples? The 7p and DIS rates are compatible in y* and at low 
Pt {y* is the 'jp cms rapidity defined with respect to the direction of the exchanged 

settings (default/DELPHI/E665): PARJ(2) = A s = (0.3/0.23/0.2), PARJ(ll) = 0.5/0.365/0.5), 
PARJ(12) = (0.6/0.41.0.6). PARJ(2) is the strangeness suppression factor, PARJ(ll) (PARJ(12)) 
is the probability that a meson containing u or d (s) quark has spin 1. 

6 For non-diffractive samples HI requires the presence of at least 500 MeV of energy in the forward 
region of the liquid argon calorimeter, 2.03 < r\ < 3.26, to remove large rapidity gap events. 
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boson). Comparison between K° rates of production in pp collisions at 200 GeV 
obtained by the UA5 collaboration! and this measurement show good agreement^. 
In contrast, the A rates in the photon fragmentation region in photoproduction are 
significantly lower than in pp interactions. This may be attributed to the presence of 
an initial baryon in pp in this hemisphere from the incoming proton or antiproton. 
A difference between photon and target fragmentation region with respect to baryon 
production has also been observed in deep inelastic [iN scattering^. 

ZEUS examines! the features of inclusive K° mesons produced with a high-i^ jet. 
The aim of the study is also to evaluate average longitudinal fragmentation function 
D(z) for K° in jets from hard photoproduction and compare it with the fragmentation 
functions! obtained in high energy quark jets in e + e~ . 



ZEUS 1994 PRELIMINARY 




Fig. 3. (a) Corrected numbers of K$ per jet as a function of jet pseudorapidity. (b) Corrected 
fragmentation function D(z). Both distributions are defined with respect to 8 GeV hadron jets. 

Events with at least one jet of E t > 8 GeV and \rf et | < 0.5 and K s with 0.5 < p t < 
4 GeV and \r]\ < 1.5 were selected. The data were corrected to take into account effects 
due to reconstruction efficiency, event selection efficiency and bin-to-bin migration. 
Comparison of the mean number of K s per jet, integrated over R < 1 (R is the jet 
radius in cone algorithm in rj and 0) are shown in Fig. [5]a as a function of rf et where 
the total number of jets was evaluated without the K s requirement. The shapes of 
the measured and Monte Carlo distributions are similar but the Monte Carlo is higher 
than data as observed also in other distributions. This discrepancy which cannot be 
explained by the size of systematic errors nor variation of MC parameters or photon 
parton densities is under investigation. 

The distributions of particles within a jet may be described by the fragmentation 
function D(z), where z is the longitudinal fraction of the jet momentum carried by 
the particle. Fig. shows D(z) for data compared to that from PYTHIA. Again 
Kg distribution is lower than Monte Carlo by 20%. The curves are upper and lower 
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limits from the phenomenological fit of Binnewies et alE The lower curve corresponds 
to gluon fragmentation and the upper curve to s and u fragmentations and covers 
the spread of the fitted D(z) values of fragmentation of u, d, s quarks using the LO 
and NLO schemes at a QCD scale of 8 GeV. The disagreement between data and 
phenomenological fit of D(z) can be due to the fact that data are corrected to the 
hadron level to avoid the strong dependence of the correction to the parton level on 
the Monte Carlo used. One should keep in mind that e + e~ data do not provide a 
sufficient constraint on g — > Kg fragmentation and in this place the ep measurement 
is important. 

2.2. Charged particles 

Latest results of the inclusive transverse momentum spectra of charged particles 
in photoproduction events in the laboratory pseudorapidity range —1.2 < rj < 1.4 
have been measured up to £><_= 8 GeV using the ZEUS detector! and in range |?7| < 1 
and 2 < p t < 12 GeV by HlM at the average jp cm. energy of W = 200 GeV. The 
inclusive transverse momentum spectra fall exponentially in the low p% region. The 
non- diffract ive data show a pronounced high p t tail departing from the exponential 
shape. Comparison to hadron-hadron collisions^ at a similar cm. energy shows 
that the photoproduction data are clearly harder and in fact are similar to pp at 
^/s = 1.8 TeV (see Fig. In Fig. |5|, the inclusive charged particle cross section da/drj 
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Fig. 4. Inclusive cross section as a function of p t 
multiplied by 10~ 4 and fitted to the dependence (1 
explains possible reasons of the difference between 



Pt spectrum is harder in 7p due to: 

• 7p is measured in 7 fragmenta- 
tion region whereas pp in the cen- 
tral region 

• 7p behaves like Vp and parton 
momenta of quarks in the vec- 
tor meson V are harder than in 
baryons 

• dominance of direct part in 7^ 
cross section at high p t can lead 
to harder scattering 

for charged particles for \ri\ < 1. pp data were 
+ ^0 -n - The text in the right part of the figure 
72? and pp distributions. 
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is shown for particles with p t > 2 GeV. A comparison is made with PYTHIA with two 
different parton distribution functions: GRV-LO(Fig. |5|a) and LACl(Fig. |5|b). None 
of the predictions match the data and they have mutually opposite slopes of the 77 
distribution. From the figure it is clear that the rapidity distribution is sensitive to the 
photon structure function — the fact which will be pursued in the final analysis of the 
data with the aim to extract the gluon structure function of the photon. The inclusive 
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Fig. 5. Inclusive cross section as a function of r\ for charged particles with pt > 2 GeV compared to 
PYTHIA with the GRV-LO(a) and LACl(b) photon parton parametrisations. 

distributions of charged particles were compared with NLO QCD calculations. Good 
agreement is observed in p t distribution in the whole measured range of pt > 2 GeV 
as well as in da/drj distribution!. 



3. Scale influence on the energy dependence of photon-proton cross section 

The energy dependence of the real photon-proton total cross section is consistently 

described at high energies by the power law a*"' oc (W 2 ) x , where W is the 7p centre 

of mass energy and the power A ~ 0.08. In contrast, the virtual photon-proton cross 

section is found to rise fasterSlii with increasing W 2 for Q 2 values larger than a few 

GeV 2 . In the HERA regime fits of the form er*»l oc (W 2 ) x lead to values of A which 
i — in 

riseat 2 ! from 0.2 to 0.4 in the Q 2 range from 1.5 to 10 3 GeV 2 (see Fig. ^]). The change 
of the slope parameter A with increasing Q 2 is associated with a transition from the 
non-perturbative 'soft' to the perturbative 'hard' regimeR The perturbative QCD 
(pQCD) can be applied in the domain where Q 2 is larger than a few GeV 2 and it 
predicts an increase of A with Q 2 . In the Reggeon Field Theory (RFI) the change of A 
can be interpreted as a reduction of screening corrections^ with increasing Q 2 . RFI 
and pQCD are two complementary approaches which successfully describe physics 
processes in different regimes - 'soft' and 'hard'. 

Whereas in DIS the hardness of the interaction (scale) can be characterized by the 
photon virtuality Q 2 , in photoproduction with Q 2 « the largest pt of a hadron in 
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the photon fragmentation region can take over the role of the scale. This investigation 
was the subject of the study disscussed in this chapter the results of which have been 
published recently^. 

In total about two million photoproduction events were selected in the energy 
range 150 < W 1P < 250 GeV. The W dependence of the cross section was measured in 
different bins of the maximal pt, ma x of the charged particle in an event. Good detector 



A 




4p 2 max /GeV 2 Q 2 /GeV 2 

Fig. 6. Slope A of the W 2 dependence of the photoproduction (a) and DIS (b) cross section as a 
function of 4p 2 and Q 2 respectively. 

acceptance and high resolution of p t measurement together with large statistics of 
photoproduction events, allow a continuous coverage of both soft and hard domains 
in the pseudorapidity range 1.1 < rj* < 3.1 (77* is measured with respect to the 
7 direction in the 7p centre of mass system). The energy dependence of the cross 
section in each pt^nax bin was fitted to the form a(W,p 2 tmax ) ~ W 2X{p ^ ] . The 
results of the fit are shown in Fig. where the slope A is given as function of ^p 2 max - 
A strong increase of A is observed with increasing p t) max which can be fitted by a 
straight line. Note that when the scale kp\ max is multiplied with a number, this 
does not change the slope of the straight line. Hence the change of A with the scale 
is significant. Comparison with the longitudinal phase space model shows that the 
phase space plays only a minor role. 

On the other hand the Pt, ma x dependence of A is described by PYTHIA. The 
model is based on the Monte Carlo version which includes both pQCD for the hard 
scattering and a Regge inspired soft interaction component. It was checked that the 
effects of different proton and photon structure function, of parton showers as well 
as multiple parton-parton scattering are small. Therefore, the change of the slope A 
is not associated with the growth of the gluon density in the proton at small x but 
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rather with the leading order matrix element of the hard process in PYTHIA. 

4. Resonances in multi-photon final state 

As mentioned already in the previous chapter the energy dependence of soft 
photon-hadron interactions at high energies are well describedtH by the Reggeon Field 
Theory provided that the exchange of Pomeron trajectory with the quantum numbers 
close to vacuum, i.e. with the intercept a(t—0) ~ 1.08 and C parity +1 is introduced. 
No direct evidence has been found at high energies for an exchange of trajectory with 
the same intercept but with C — — 1 called Odderon trajectory. At HERA the exis- 
tence of an Odderon trajectory could result in the production of pseudoscalar mesons 
with C = +1 in addition to their production via the two photon fusion process. 
Schafer et al.0 suggest to measure the cross section ratio for exclusive production 
of 77's and r) n s since the sensitivity to the Odderon contribution differs for 7] and r)'. 
Assuming the coupling to the meson for photon-Odderon fusion being proportional 
to the product of the electric charge and the baryon number of the quarks composing 
the meson, 77 - being predominantly in the octet state with respect to flavour SU (3) - 
gets larger contribution compared to rf being dominantly the singlet. The data were 
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Fig. 7. Invariant mass spectra of 27's: (a) 7r° mass region for photon energies E 7 > 1 GeV after 
background subtraction (b) 77 mass region for photon energies E 1 > 1.5 GeV after background 
subtraction (c) events with a vertex and E 1 > 0.7 GeV (d) events without a vertex and E 1 > 
0.7 GeV. 

taken in 1996 with the HI detector and amount to an integrated luminosity of about 
5 pb _1 . The lead/scintillating fibre calorimeter SPACALO has a sampling term of 
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7%/yE/GeV and a position resolution of about 5 mm at 1 GeV. Together with the 
threshold « 100 MeV this allows for a precise measurement of mesons decaying into 
multi-photon final states, e.g. of 7r° and rj mesons in their 2 / y decay mode, as shown in 
Fig [7[ The central values of the Gaussian fit amount to 136 MeV for 7r° and 535 MeV 
for rj. The measured width of the invariant mass distribution is a « 13 MeV for 7r° 
and a ~ 31 MeV for t], demonstrating the good performance of the detector. 

The acceptance of the positron calorimeter of the luminosity system restricts the 
photon energy in the laboratory system to lie between 8 and 20 GeV and to a photon 
virtuality Q 2 < 0.01 GeV 2 . The selected events consist of the scattered positron, 
two photons within the angular acceptance of the SPACAL (—3.8 < rp < —1.4, rf 1 
is photon laboratory rapidity), and an elastically scattered proton which escapes 
undetected in the main detector. 

Selected events with a reconstructed vertex, i.e. with at least one charged track in 
the central detector, result in the invariant mass spectrum of photon pairs as shown 
in 0c. A clear signal of 7T° and rj mesons are observed, no indication for a signal of 
the rf meson at the mass 958 MeV is found. By requiring no charged particles in the 
central detector (events without a reconstructed vertex) is compatible with exclusive 
meson production. It leads to the invariant mass spectrum of the photon pairs as 
shown in [7]d. No signal of exclusively produced rj and rf mesons is observed nor of 7r°. 
The resonance-like structure at an invariant mass of about 800 MeV is interpreted as 
being due to elastic photoproduction of u mesons and their subsequent decays into 
the 77T° final state, where only two out of three photons are detected. 

The non-observation of rj mesons is compatible with the expectation from a pure 
two photon fusion process, given the small acceptance of the detector. This excludes 
a large contribution from Odderon induced process with cross section substantially 
larger than the two photon fusion cross section. In the case of rf the small branching 
ratio of 2.12% further reduces the expected number of events by an order of mag- 
nitude. The analysis continues to quantify the results into the cross section upper 
bounds for exclusive production of 7r°, rj and rf mesons. 

5. Observation of isolated high E t photons at HERA. 

The ZEUS detector identified a class of events showing the characteristics of high 
E t isolated photons termed here as "prompt" to distinguish them from those produced 
via particle decays. This is the first observation of prompt photons at 7p centre of 
mass energies an order of magnitude higher than those previously employed. In the 
kinematic range covered by HERA detectors, the prompt photons are produced pre- 
dominantly in two hard subprocesses0 which in LO are represented as the Compton 
scattering ^yq — > ^q in the direct channel and by qg — > ^q in the resolved channel. 
NLO QCD adds additional processes where hard photon is radiated from outgoing 
quarks. Due to the electromagnetic nature of the photon the prompt photon pro- 
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duction is suppressed approximately by a factor a with respect to hard processes 
involving partons only which represent the main source of background. Therefore, a 
sophisticated data analysis must be done to isolate events with prompt photons from 
background. A particular virtue of prompt photon processes is that the observed 
final-state photon emerges from the QCD process directly, without the intermediate 
hadronisation of final state partons into jets. 

The data correspond to an integrated luminosity of 6.36 pb^ 1 collected by the 
ZEUS detector in 1995 in e + p collisions with energy E e = 27.5 GeV with protons 
of energy E p = 820 GeV. Photons were detected in the barrel part of the uranium- 
scintillator calorimeter in its electromagnetic section (BEMC) which restricts the 
photon pseudorapidity range to —0.75 < rf< < 1.0. A standard ZEUS electron finding 
algorithm was used to identify photon signal with 5 < El < 10 GeV. A cone 
jet finding algorithm was used to identify jets with E\ et > 5 GeV, pseudorapidity 
— 1.5 < if et < 1.7 and unit cone radius, using entire ZEUS calorimeter system. Events 
with an observed DIS electron were removed, restricting incoming photon virtualities 
to Q 2 < 1 GeV 2 . Isolation cone criteria were applied on the photon candidate to 
reduce possible background from radiation of electrons or positrons and from energy 
fluctuation of jets into hard n° or rj. Two shape-dependent quantities were studied 
in order to distinguish photon, 7r° and 77 signals. These were the energy weighted 
mean width of the BEMC cluster in z direction z W idth an d the energy fraction of the 
most energetic cell in the cluster F max (Fig. |S]a) of the photon candidate. The relative 
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Fig. 8. (a) Distribution of F max for prompt photon candidates and 7, 7r°, 77 fitted Monte Carlo. z W idth 
given in BEMC cell units (b) x™ eas of prompt photon after background subtraction. (•) data points, 
histogram - Monte Carlo: (■ • ■) quark radiation, ( ) q rad+resolved, ( — ) q rad+resolved+direct. 

amounts of 7, ir° and rj from single particle Monte Carlo distributions were obtained 
from fits of z width and F max to data. To get distribution of physical quantities, a 
background subtraction is applied on the assumption that the data may be expressed 
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as a sum of the photon signal plus neutral meson background. The result for the 
fraction x 1 of the incoming photon momentum that contributes to the production of 
the high Et photon and jet x™ eas = X)j e t+7 event 

{E-p z ) is shown in Fig. gb. 

The Monte Carlo distributions are normalised to the same integrated luminosity as 
data. A strong peak near x™ eas = 1 corresponds to the Compton process. It is not 
possible to make conclusions concerning the presence of a resolved photon component. 
The background distribution is consistent with zero for x™ eas > 0.9. The quoted 
cross section for the process ep — ► e + jprom.pt + J e ^ + X with x™ eas > 0.8 is 
15.3 ± 3.8(stat) ± 1.8(syst) pb in the kinematic range defined above. The systematic 
error is dominated by 8% contribution due to the uncertainty in the calorimeter 
calibration. The cross section is in an agreement with the NLO QCD calculation^. 

6. Conclusion 

Inclusive particle production provides a complementary view on hard scattering 
to the commonly more used jet production. Selecting events with a hadron with p t 
at least 1 — 2 GeV, the inclusive spectra are calculable in QCD and are sensitive to 
structure functions of beam projectiles — photon (real or virtual) and proton. It 
will be the task of the coming years of HERA hadron physics to use inclusive spectra 
of hard hadrons to obtain the gluon structure function of the real photon. Events 
with prompt photons have a sensitivity to the quark content of the photon. With the 
increase of HERA luminosity they will supplement the inclusive spectra of charged 
particles in the study of partonic structure of the photon. 
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